Vortex Pinning and Non-Hermitian Quantum Mechanics 
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A delocalization phenomenon is studied in a class of non-Hermitian random quantum-mechanical 
problems. Delocalization arises in response to a sufficiently large constant imaginary vector poten- 
tial. The transition is related to depinning of flux lines from extended defects in type-II supercon- 
ductors subject to a tilted external magnetic field. The physical meaning of the complex eigenvalues 
and currents of the non-Hermitian system is elucidated in terms of properties of tilted vortex lines. 
The singular behavior of the penetration length describing stretched exponential screening of a per- 
pendicular magnetic field (transverse Meissner effect) , the surface transverse magnetization, and the 
trapping length are determined near the flux-line depinning point. 
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I. INTRODUCTION 

Although Hamiltonians must be Hermitian in conven- 
tional quantum mechanics, non-Hermitian operators do 
appear in other physical contexts: the time evolution 
of non-Hermitian Liouville operators can describe vari- 
ous nonequilibrium processes; the transfer matrix 
of two-dimensional asymmetric vertex models leads to 
non-Hermitian Hamiltonians for quantum spin chains. ^] 

In the present paper, we investigate localization phe- 
nomena in an especially simple class of random non- 
Hermitian Hamiltonians. Although non-Hermitian, our 
problem is sufficiently close to conventional quantum me- 
chanics that it will be convenient to use a quantum lan- 
guage to describe the results. Specifically, we show that 
a delocalization transition occurs (even in one and two 
dimensions) in the following one-body Hamiltonians in d 
dimensions: first, the Hamiltonian in continuum space. 



n 



jp + ig) 

2m 



V{x), 



(1.1) 



where p = (h/i)d/dx is the momentum operator and 
V{x) is a random potential; second, the second-quantized 
lattice Hamiltonian, namely the non-Hermitian Anderson 
model on a hypercubic lattice, 

+ (1-2) 

X 

where the vectors {e^} are the unit lattice vectors, the 
{bx,bx} are (boson) creation and annihilation operators. 



Vx is a random potential. In both of the Hamiltonians, g 
is a non-Hermitian external field. Although the bulk of 
our discussion will concentrate on the pro perties of the 
single-particle Hamiltonian (1.1) and (1.2), many of our 
results will be relevant for interacting many-body boson 
problems, provided that we forbid double occupancy of 
eigenstates in the localized regime. Interaction effects 
in both the localized and delocalized phases will be dis- 



cussed in Sec. VIII 



We can regard the non-Hermitian field as an imaginary 
vector potential. Models with a real gauge field A would 
be written in the Hermitian forms 



n' 



2m 



Vix) 



(1.3) 



and 



+Y.^xbibx. (1.4) 



In two dimensions with spatially varying A — A{x), 
these Hamiltonians describe the quantum Hall system, 
where some of the localized states of the A = G case 
are delocalized in the presence of the gauge field. ^ We 
obtain the non-Hermitian Hamiltonians ( |1.1| ) and (1.2) 
by replacing —eA{x) with a constant, ig. In this non- 
Hermitian case, we show that all eigenstates can be de- 
localized (even in one dimension) for large g. 
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FIG. 1^ A situation described by the classical Hamilto- 
nian (1.7) in (1 + 1) dimensions. The wavy line indicates 



a tilted string with single-valued trajectory x{t) subject 
to thermal fluctuations. 
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FIG. 2. Flux-line depinning from a columnar defect 
in a superconductor: (a) the bulk part of the flux line 
is pinned by the defect, and the transverse magnetiza- 
tion has contribution only from the surface transverse 
magnetization M^g] (b) the flux line is depinned, and 
the "surface" transverse magnetization now has a macro- 
scopic value Mj^s cx L-^^ where is the thickness of the 
sample. 



This problem has direct physical relevance when we 
map the non-Hermitian quantum system in d dimensions 
to a problem of classical equilibrium statistical mechan- 
ics in (d -t- 1) dimensions. Consider a matrix element of 
the imaginary-time-evolution operator, 



(1.5) 



where the Hamiltonian Ti. is given by (1.1). The bra and 
ket vectors are the initial and final states, respectively. 
Using the standard path-integral scheme, we can express 
the above matrix element as the partition function 



(1.6) 



where J Dx denotes summation over all possible world- 
line configurations £c(t), and Ed is the energy of a clas- 
sical elastic string in (d + 1) dimensions, given by 
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(1.7) 



The energy Ed is of course the imaginary-time action for 
the equivalent quantum problem. The Planck parame- 
ter h is interpreted as the temperature of the classical 
system. 

The above classical system expresses the following 
physics; see Fig. |l|. Consider a thermally fluctuating 
string put on a r ando m washboard potential. The first 
term of the Eq. ( |l.7[ ) is optimized when the string is 
tilted from the r direction by the angle = tSLn~^ (g/m). 



On the other hand, the potential tries to make the string 
parallel to the r direction, to take advantage of partic- 
ularly deep valleys in the random potential. We might 
thus expect a depinning transition of the string from the 
random washboard potential with increasing g. 

We suspect that there are many physical realizations 
of the above statistical-mechanical problem. For an ap- 
plication to population biology, see Ref. [t) . The present 
study is motivated particularly by a depinning transition 
of magnetic flux lines in type-II superconductors with 
quenched disorder. It has been recognized that impuri- 
ties and defects play an essential role in the mixed phase 
of type-II superconductors. Flux lines, subject to elec- 
tromagnetic forces due to applied currents, can move and 
dissipate energy, thus destroying the superconductivity, 
unless they are pinned by impurities or defects. In par- 
ticular, extended defects such as columnar defects and 
twin boundaries pin flux lines parallel to them efficiently, 
and thereby expand the region of zero resistance dramat- 
ically; see Refs. @- 0. 

The situation described by the classical Hamilto- 
nian (1.7) emerges when the external magnetic field 



generating the flux lines is not parallel to the extended 
defects; see Fig. ^. We may expect that the bulk part of 
the flux line remains pinned by the extended defects when 
the transverse component of the magnetic field, H±, is 
weak. This phenomenon is called transverse Meissner ef- 
fect; the system exhibits perfect bulk diamagnetism 
against H± for a slightly tilted magnetic field. When 
we increase H±, i.e. tilt the magnetic field further rela- 
tive to the extended defects, a depinning transition takes 
place at a certain strength of the transverse field, namely 
H±c. We note that randomness in the potential is essen- 
tial for the transverse Meissner effect to work at finite 
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FIG. 3. A schematic phase diagram of high-temperature 
superconductors with columnar defects. |]10[ The abscissa 
indicates the transverse component of the external mag- 
netic field. 



Anderson transition and Mott's variable-range hopping. 
Our problem is perhaps the simplest non-Hermitian gen- 
eralization of these phenomena. The inverse localization 
length of a wave function, when it delocalizes, is pro- 
portional to the critical transverse field H±c- As H ± 
increases, a precursor of the flux-line depinning appears 
first near the surface of the superconductor. Slightly be- 
low the critical point, the ends of the pinned flux line be- 
gin to tear free from the pinning center. Thus the trans- 
verse Meissner effect first breaks down near the surfaces, 
and the typical depth of the region where this bending 
takes place is a "penetration depth" [different from the 
penetration depth associated with the underlying super- 
conducting material) which accompanies the transverse 
Meissner effect. A Mott variable-range- hopping descrip- 
tion of this depinning will lead us to the conclusion that 
the penetration depth r* diverges at the transition point 
Hi^c as 



for the random columnar defects, and 



(1.8) 



temperatures. For a periodic potential, thermally acti- 
vated kinks always lead to extended eigenf unctions in the 
equivalent quantum problem (Bloch's theorem!), and the 
effect vanishes. See Ref. for a discussion. 

Nelson and Vinokur [|10| have described flux-line pin- 
ning by means of the path-integral scheme. They mapped 
a phenomenological Hamiltonian for flux lines with ex- 
tended defects onto interacting bosons in 2 -f 1 dimen- 
sions with point impurities. A low-temperature phase 
with flux lines localized on extended defects was related 
to the Bose-glass phase ||l2| of disordered boson sys- 
tems. When the applied magnetic field is tilted away 
from the extended defects, the corresponding boson prob- 
lem acquires the non-Hermitian field g. |^,|^ The flux- 
line depinning transition due to the magnetic-field tilt 
corresponds to a delocalization transition of the non- 
Hermitian quantum problem. The expected T-H± phase 
diagram of the flux-line system (with J?||, the field com- 
ponent parallel to the columns, held fixed) is shown in 
Fig. ^. The transition to a flux-liquid or "superfluid" 
phase with increasing temperature for H ±^ = is of the 
usual Bose-glass type. [|2| In the present paper, we shall 
focus on the transition curves for H ± ^ 0. Hence we fix 
the temperature and follow the trajectory of increasing 
H±, indicated in Fig. || as a thick arrow. We first de- 
scribe the delocalization transition of the non-Hermitian 
quantum problem, and next translate the results into the 
language of flux-line depinning. We mainly discuss the 
one-body (or one-flux-line) problem of the non-Hermitian 
quantum mechanics. Some of our results were summa- 
rized in a previous publication, js] 

One advantage of discussing the flux- line depinning in 
terms of an equivalent "quantum" system comes from 
the fact that many ideas are already available concern- 
ing localization in Hermitian random systems, such as the 



(1.9) 



for the random twin boundaries. The pinning becomes 
ineffective above the critical field. The direction of the 
bulk transverse magnetization jumps at the critical point, 
while the surface transverse magnetization diverges in 
both cases according to 



M±s ^ {H±c — H\ 



(1.10) 



We find stretched- exponential relaxation of the perpen- 
dicular magnetization to zero as one goes deep into the 
bulk of the sample from the surface. 

The penetration depth r* may be observed experimen- 
tally by measuring the response to an AC magnetic field 
superimposed on the DC field H ^\ only the dangling 
ends of the flux lines will respond to the AC field, and 
the length of these line segments diverges at the transi- 
tion. It would be interesting to check our prediction of 
a diverging penetration depth near the Bose-glass tran- 
sition with such an experiment. 

A preliminary account of this work appeared in Ref. . 
For a related problem which arises from the physics of 
charge density waves, with some results applicable to vor- 
tex lines, see Ref. [Q. See Ref. |15 for a related problem 
in fluid mechanics. 

The plan of this paper is as follows. Section || de- 
scribes basic relations between the flux-line syst em and 
non-Hermitian quantum mechanics. In Sec. |lll| , we ex- 
plain how the delocalization transition occurs in the non- 
Hermitian random system. We present solutions of non- 
Hermitian systems in Sees. |rv|-|Vl[ The prediction s for 
the transverse Meissner effect are derived in Sec. VH. 



Section VIII discusses the effect of interactions 
delocalized phase in both 1 + 1 and 2 



in the 
1 dimensions. 
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II. FLUX-LINES AND NON-HERMITIAN 
QUANTUM MECHANICS 

In this section, we review the basic correspondence 
between vortex hnes in superconductors and the non- 
Hermitian quantum system (1.1). We then describe what 
the current and the wave functions of the non-Hermitian 
system mean for vortex trajectories. 



A. Path-integral mapping 

We start with the energy of a flux hnc in {d + 1) 
dimensional type-11 superconductor: |lO| 




dx 



dr 



Vix) 



Air 



(2.1) 

This is a phenomenological Hamiltonian for the flux-hne 
system where all columnar defects or twin boundaries are 
parallel. Here r denotes the coordinate parallel to the 
extended defects, while x denotes the d-dimensional co- 
ordinates perpendicular to the defects. We describe the 
flux line as a single- valued function x{t) by neglecting 
overhangs of the flux line. The thickness of the super- 
conductor in the r direction is denoted by L^- 

The first kinetic-energy-like term of the integrand 
in ( ^.l[ ) comes from the harmonic approximation to the 
energy increase due to local tilt of the flux line. The 
"mass," or tilt modulus is denoted by ii. When us- 
ing quantum language to describe this problem, we shall 
make the correspondence 

ei < — > m. (2.2) 

The potential V is generated by columnar defects or 
twin boundaries. Each columnar defect is specified by a 
point Xk in the x plane. We assume that the columnar 
defects are distributed randomly: 



M 

V{x) ==^l/i(a;- 

k=l 



(2.3) 



where M denotes the number of the columnar defects, 
Vi is the potential of an individual defect, and Xk is the 
random position of the fcth columnar defect. In the case 
of twin boundaries, we consider only the situation where 
all twins are parallel to each other. By projecting out the 
degree of freedom perpendicular to r and parallel to the 
twin boundaries, |jl^,|l6| we can reduce the dimensional- 
ity of the a;-space by one. (We neglect effects due to an 
enhanced concentration of point disorder which may oc- 
cur in twins.) Thus each twin boundary acts like a line 
defect in 1 -I- 1 dimensions. In short, we can set d = 2 for 
columnar defects and d = 1 for twin boundaries. 

is due to the 



The third term of the integrand in ( 2.1 



transverse component of the magnetic field H± . The fiux 
quantum 27r?ic/(2e) is denoted by (j)Q. In the absence of 




FIG. 4. One flux line (wavy curve) induced by the field 
and interacting with columnar defects in a cylindrical 
superconducting shell. The transverse component H± is 
generated by the current / threading the ring. 

a pinning potential, this term is optimized energetically 
when the flux line tilts to follow the external fleld. 

Figure ^ shows a flux-line system described by the 
above phenomenology in (1-1-1) dimensions with colum- 
nar pinning. We assume periodic boundary conditions 
in the x direction, although these could be replaced by 
reservoirs of flux lines at the edges of the sample. In this 
geometry, the transverse component H ± is generated by 
the current threading the cylinder of a thin superconduc- 
tor. As we shall see, the physics resembles the Aharonov- 
Bohm experiment carried out on a mesoscopic quantum 
ring with, however, an imaginary "flux" through t he ri ng. 

The partition function of the flux- line system (2.1) is 
given by 



Z = Vx exp — 



(2.4) 



where J T>x denotes integration over all possible config- 
urations 3; (r) of the flux line. Upon identifying Eq. (2.4) 
with Eq. ( |l.6[ ), we are led to a correspondence between 
the flux-line syst em ( ^.l| ) and the non-Hermitian quan- 
tum system ( [l.l[ ) summarized in Table |. In particular, 
the thermal fluctuations of the flux-line system are equiv- 
alent to zero-point motion of the quantum system, and 



kBT 



(2.5) 



The tilt term in (2T) is related to the non-Hermitian fleld 
in (0) by 

9^^- (2.6) 
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TABLE I. The correspondence between the (d + l)-dimensional flux-line system and the d-dimensional quantum system. 



Partition function ( |2.4[ ) Partition function (|l^ 

Thickness Lr Imaginary time r 

Temperature ksT Planck parameter h 

Tilt modulus ii Mass m 

Tilting fleld (f>oH±/{4:n) Non-hermitian fleld g 

Tilt slope dx/dr (oc Transverse magnetization) Im J, where J = current 

Probability distribution at r = (L^ — > oo) Left-eigenvector i}^{x) 

at r = Lt Right-eigenvector ip^lx) 

at r = Lr/2 tp^^i^ 



Boundary conditions at the bottom and top surfaces of 
the superconductor are represented by the initial and fi- 
nal vectors of Eq. (1.5) in the forms 



dx (x\ 



and 



1^0 



dx \x) 



(2.7) 



These boundary conditions in the imaginary-time direc- 
tion mean that we integrate freely over the points where 

vortex lines enter and exit the sample. 

Although the derivation of the quantum proble m (|l.l| ) 
from the classical elastic string described by Eq. ( |2.l| ) is 
straightforward in the path-integral scheme, it is also use- 
ful to think in terms of Galilean b oost s. p^ , p^ We can 
eliminate the additional term in ( ^.l[ ) by the Galilean 
transformation 



x' = X — vt, 
t' = t. 



where the velocity of the moving frame is 



g_ 

im 



(2. 



(2.9) 



Note that the velocity is imaginary, because of the imag- 
inary time t = ir in the equivalent quantum problem. 
The above Galilean transformation changes the kinetic 
term of the Hermitian Hamiltonian to the form 



2m 



2m 



(2.10) 



which results in the kinetic term in (1.1). 

We do not explicitly include interactions between flux 
lines in t he H amilt onia n (2.1) nor in the quantum Hamil- 
tonians (IT) and (L2). However, strong short-range re- 
pulsive potential may be treated approximately by for- 
bidding multiple occupancy of localized levels of the 
quantum Hamiltonian. Since the particle mapped from 
the flux line follows the Bose statistics, |l^ we may treat 
many-flux-line problems as hard-core bosons in a random 
potential. For d = 1, hard-core bosons are equivalent to 
fermions, and hence we can use the Pauli principle di- 
rectly to treat many-body problems in this case. See 
Sec. VH] for a discussion of interaction effects using a 



boson formalism in d dimensions. 



B. Imaginary current and the transverse 
magnetization 



We show in the following that the imaginary part of the 
current of the non-Hermitian system yields the average 
slope of a vortex trajectory as it crosses the sample. 

Since the non-Hermitian field g acts like an imaginary 
vector potential, we define the current operator as 



J 



. dH p + ig 
dg m 



where 



' para 



P. 

m 



and 



' para 



J dia 



(2.11) 



(2.12) 



Note first that the expectation value of the position of 
the flux line at the distance r from the bottom surface 
of the superconductor is given by 



1 



(2.13) 



where the partition function Z is defined by (1.5) and a;°P 



is the position operator acting on position eigenstates \x) 
such that a;°P|a;) = x\x). Note also the useful commuta- 
tion relation 



[H,a;°P] = -ihJ, 
which leads immediately to 



(2.14) 



^ ~i (J)^ = Im (J) 



(2.15) 



In the last line of Eq. ( 2.15 ), we used the fact that the 
expectation value of the current is pure imaginary in the 
present problem. The current carried by an individual 
eigenstate is not necessarily pure imaginary, where the 
current of a state e„ is defined by 



^ _.den 

dg 



(2.16) 
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However, a purely imaginary expectation value of the 
current arises owing to a pairing property of eigenval- 
ues and eigenfunctions: if there is a complex eigenvalue 
e with a right-eigenvector ?/;^ for this non-Hermitian 
problem, there is also the eigenvalue g* wi th the right- 
eigenvector (V'^)*- When we sum Eq. ( 2.16| ) over all the 
eigenstates to calculate the expectation value, the real 
part of the current of a state cancels out in pairwise fash- 
ion. 

The transverse component of the magnetization per 
vortex line at depth t below the surface of the sample is 
proportional to the tilt slope. 



d 



(2.17) 



The total transverse magnetization is proportional to the 
net displacement of the flux line per unit length between 
the bottom surface (t = 0) and the top surface (t = Lr): 



M_L = nyy^f mj_{T)dT = ^^Im / dr (J)^ 
(poriv 



'0/ ' 



(2.18) 



where n-u is the density of vortices. As we show below, 
the last factor in the brackets is given by 



dg 



(2.19) 



Equation (2.18) is an indicator of the dclocalization tran- 
sition; the transverse magnetization must remain finite in 
the limit L,- ^ oo in a depinned phase, while it must van- 
ish macroscopically in a pinned glassy phase, which is the 
transverse Meissner effect. In fact, we show later that the 
vortex tilt m±{T) appears only close to surfaces in the 
pinning regime, and th us M ± vanishes as Lr ^ oo. In 
the limit Lr ^ oo, Eq. (2.19) becomes 



(2.20) 



where Egg denotes the ground-state energy of H(g). (The 
ground state is defined by the lowest real part of the 
eigenenergy. |18|) Hence, in the pinning regime, the 
ground-state energy must be independent of g. In the 
depinning regime, —dsgs/dg is nonzero and equals the 
mean tilt slope of th e corresponding flux line. 

To prove Eq. ( 2.19 ), we proceed as follows. As is stated 
above, the non-Hermitian field g acts like an imaginary 
vector potential. Because of this, we can gauge away 
the non-Hermitian field by applying the imaginary gauge 
transformation to the non-Hermitian Hamiltonian (|l . ll) : 



u~'nig)u^nio), 



where 



U = exp 



gx°P 



(2.21) 



(2.22) 



(This is lust another expression of the Galilean transfor- 
mation ( 2.1C| ).) Therefore, we have 



(2.23) 



and 



d 



dg 



^-Lrn{g}/h^op^ 



(2.24) 



Upon taking thermodynamic averages (ex pecta tion val- 
ues in quantum language), we obtain Eq. ( ^.19 ). 

A conductivity-like transport coefficient may be de- 
fined by 



d_ 



The quantity defined for each eigenstate 



t(») = 



(2.25) 



(2.26) 



is a measure of the stiffness of the eigenfunction. In 
many-body problems, a^i, is related to the superfluid 
density of the relevant boson system. p9l 



C. Probability distribution of a flux line and 
eigenvectors of the non-Hermitian Hamiltonian 

We now outline the correspondence between the prob- 
ability distribution of a flux line and the eigenvectors of 
the quantum system. The flux line fluctuates because of 
thermal excitations. This thermal agitation is described 
by quantum fluctuations of the corresponding fictitious 
non-Hermitian quantum particle. The probability distri- 
bution of the flux line due to thermal fluctuation is thus 
related to wave functions of the quantum system. 

The probability distribution of the flux line at the dis- 
tance r from the bottom surface is given by |l0| 

P(a;; r) = - (V'/j e-(^^-^)«/'' \x) {x\ e'^^/'^ , 



so that we can rewrite Eq. ( 2.13| ) as 



{x)^ = J xP{x;T)d'^x. 



(2.27) 



(2.28) 



Let us calculate P{x] r) in the limit Lt — > oo. In this 
limit the exponential operator exp(— L^7Y/?i) can be ap- 
proximated by 

hm e-^-^/"- |^g,)e-^-'-/'"MV'g«l- (2-29) 
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Note here that left and right eigenvectors may be ex- 
pressed as 

^^{x) = (cc I V) and i^^ix) = (tA | a;) . (2.30) 



Upon combining Eq. ( |2.27|) and Eq. (|2.29| ), we find 



and 



^Hm^P(.;0) = -^^^^^, 

ijg{x) 

Urn P(x;Lr)^-r--si—m—i 



_hm P{x;Lr/2)=i;^^{x)^^,{x) 



(2.31) 
(2.32) 

(2.33) 



Since the ground state wave function is positive definite, 
the distributions (2.31) and ( 2.32| ) are weU-defined. 

We remark on the nature of the hmit in Eq. ( [2.29 ). 
For one-impurity syste ms (se e Sec. |^), corrections to the 
right-hand side of Eq. ( 2.29D are of order exp(— Lt-Ae), 
where Ae is the energy gap above the ground state. For 
random systems, the argument in Sec. VI] yields correc- 



tions which vanish as a stretched exponential function of 

Lr. 

For a finite density of fiux lines, we can treat a strong 
short-range repulsive potential approximately by forbid- 
ding multiple occupancy of localized levels of the quan- 
tum Hamiltonian. [ 1^ We then fill up the localized levels 
in order of increasing energy up to a certain level e = fi. 
The system is now characterized by an average chemical 
potential /i — /i(i?||), where is the external magnetic 
field along the r axis. This chemical potential controls 
the flux-line density and separates occupied and unoc- 
cupied levels. Under this assumption, we can estimate 
the probability distribution of the most weakly bound 
flux line by using the state at the chemical potential, ^p^, 
instead of V'gs in the above analysis. 



III. DELOCALIZATION TRANSITION IN 
NON-HERMITIAN QUANTUM MECHANICS 

We now sketch the mechanism of the delocalization 
in the non-Hermitian quantum system described above. 
Consider first localized states in a small transverse field 
g. Assume that we know the eigenvalues e„ and the 
eigenf unctions tpn{x) for g = 0. We can use the imagi- 
nary gauge transformation ( 2.22| ) to determine the eigen- 
functions for small g. The right-eigenfunctions and the 
left-eigenfunctions are given by 



i^n{x;g) 

i'nix-g) 



Uipn{x;g ^ 0), 



(3.1) 
(3.2) 




FIG. 5. The wave function ( p.5| ) or ( [4.7|) , the ground 
state of a single impurity with < g < in one dimen- 
sion. 



U = exp 



gx^ 



and U 



exp 



gx'^ 



(3.3) 



The corresponding energy eigenvalue e„ is then un- 
changed provided that this imaginary gauge transforma- 
tion is applicable. However, the above wave functions ip^ 
and t/)^ may diverge as |a;| — > cx) and become unnormal- 
izable. [ pT[ If we assume an asymptotic form of the wave 
function for g = of the form 



ipn{x;g = 0) 



exp {-Kn\x - Xr- 



(3.4) 



the condition for non-divergence of -0^ and ip^ in \x\ — > 
cx) is Igl < hKn- Then the normalized wave functions are 
approximately given by (see Fig. H) 




r(d)r!rf 



exp 



exp 



g {x- Xn) 



g {x- x„) 



Kjy). X Xr 



(3.5) 



(3.6) 



where a;„ is the localization center for g — 0, Qd the 
total solid angle of the d-dimensional space, and the nor- 
malization condition is 



^^V'^rf'^a; = 1. 



(3.7) 



where 



We naturally regard the point \g\ the delo- 

calization point of the nth state. Since the eigenfunctions 
are proportional to the probability distribution of a flux 
line a t the surfaces of the superconductor (see Eqs. ( 2.31 ) 
and ( 2.32 )), we interpret the divergence in |a;| — > cxd as 
dcpinning of the fiux line. Because localized states with 
finite localization length are stable against a small g field, 
flux-line pinning is robust against small transverse field: 
this is the transverse Meissner effect. [|l0|| 

The diverging tilt slope in Eq. ( 1.10| ) follows, be- 
cause the localization length of the right-eigenfunction is 
(k„ — \g\/h)~^ in the direction of g, and (k„ -|- \g\/K)~^ 
in the opposite direction. The former diverges as \g\ 
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hKn , hence H± c = 4 7r?tKn / t^p in Eq. ( 1.1C| ). According 



to the relation ( 2.18 ), the surface displacement or "lo- 
calization length" of the flux line is proportional to the 
surface transverse magnetization AI±s. 

We need to specify boundary conditions in order to ob- 
tain a well-defined wave function in the thermodynamic 
limit. Consider for simplicity the periodic boundary con- 
ditions 



^^{L,/2,y,...)=^^i~L,/2,y,...), 



(3.8) 



where we put the x axis parallel to g. (The one- 
dimensional periodic system is realized in the setup 
shown in Fig. ^.) A wave function of the form ( |3.5| ) 
localized at the origin has a mismatch of the factor 
exp[— (k„— |g|/?i)La;] at the boundaries x = ±La;/2. This 
mismatch is exponentially small in the pinning regime 
\g\ < hKn- Hence the wave functions ( |3.5| ) and 
are excellent approximations when we take the thcrmo- 
dyna mic limit, imposing the periodic boundary condi- 
tions (3.8). 

In the depinning regime \g\ > TiKn, however, we need to 
change the wave functions drastically to meet the bound- 
ary conditions. To see what happens in this regime, con- 
sider the limit in which the random potential V{x) may 
be neglected: 



(p + ig) 

2m 



(3.9) 



(We argue later that the potential in fact has only per- 
turbative effects for l arge g when d = 1.) The periodic 



boundary conditions (3 
extended Bloch wave 



are then satisfied only by the 



^^{x) = eyip{ik ■ x) 



(3.10) 



with — 2n^TT / L^^ where is an integer and L^, is the 
system size in the Xi, direction. The eigenvalue of this 
wave function is complex: 



e{k) 



{hk + ig) 
2m 



(3.11) 



We can interpret the appearance of the imaginary part 
of the energy in the following way: A depinned flux line 
in a periodic system has a spiral trajectory, and hence 
periodicity in the imaginary time direction. (See Fig. |^.) 
Because the single-line partition function at "time" r 
may be written Z{t) = X)™ "^"^n (^)^~^"^^''' where the 
constants c„ are coefficients depending on "initial con- 
ditions" at, say, the bottom surface of the sample, the 
period of motion in the imaginary-time direction associ- 
ated with the nth eigenstate is given by ?i/Im £„, where 
n denotes the wave function describing the depinned flux 
line. Thus a complex energy, as well as an imaginary 
part of the current ( 2.19| ), is an indicator of the depin- 
ning transition in periodic systems. 




FIG. 6. A spiral trajectory of a flux line in the system of 
Fig. |4| reflecting the imaginary eigenvalues, which appear 
in non-Hermitian systems when boundary conditions al- 
low nonzero tilt of vortex lines across the system. 



We use periodic boundary conditions in all directions 
perpendicular to the imaginary-time axis throughout this 
paper. The effect of "free" boundary cond itions in the 
time-like direction is discussed in Sec. 



VII. For super- 



conducting flux lines, unusual geometries such as Fig. ^ 
would be required to exactly implement periodic bound- 
ary conditions in the direction parallel to the tilt field 
g. For conventional slab-like superconducting samples, 
tilt causes extra vortices to enter and exit the sample on 
the sides perpendicular to g. These boundary conditions 
correspond to reservoirs and sinks of particles at oppo- 
site ends of the sample in the quantum analogy. Because 
periodic boundary conditions mimic the effects of such 
sources and sinks, we expect that the basic results of this 
paper apply to slab-like geometries as well. 



IV. ONE-DIMENSIONAL ONE-IMPURITY 
PROBLEM 



In the present section, we analyze the non-Hermitian 
Hamiltonian (1.1) analytically in one dimension with a 
single impurity. The calculation is instructive as well as 
of physical relevance for vortex lines in the case of a twin 
boundary. Throughout this section, we assume without 
loss of generality that the non-Hermitian field is applied 
in the positive x direction, i.e. g > 0. 
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A. Exact solution: a point impurity 

To determine the effect of an isolated impurity, we cal- 
culate the right-eigenvectors of the Schrodinger equation 
T-lil}^{x) = eip^{x) with the Hamiltonian defined by 



n 



{p + 

2m 



VoSix) 



2m 



Vo5{x) (4.1) 



with Vo > and the periodic boundary conditions, 



(4.2) 



The details of the solution are given in Appendix ^. 
There is a unique localized ground state for g < gc, where 
the critical field is 



mVo 



(4.3) 



(Recall that the ground state is the state with the lowest 
real part of the eigenenergy.) This solution is the only lo- 
calized state in the entire spectrum. There are extended 
excited states for g < gc as well as extended ground and 
excited states for g > gc- For later use in the next subsec- 
tion, we write down the solutions of the extended states 
including finite-size corrections. 

First, the ground state for g < gc has the energy 



2m 



mV^ 
' 2f? 



with 



mVo 



(4.4) 



(4.5) 



Note that the ground-state energy is independent of g 
and equal, in particular, to the ground-state energy in 
the Hermitian case g = 0. This observation is consistent 



with the pinning criterion given below Eq. (2.20). The 
current carried by the ground state therefore vanishes; 



Im Jjj, 



dg 



0. 



(4.6) 



The corresponding eigenvector is a localized function 
written in the form 




FIG. 7. The g dependence of the ground-state energy 
(solid line) and the first-excited-state energy (dashed 
line) of the one-dimensional one-impurity problem. 



with 



/tox = I + ^ In + 0(L-'), 

n gc - g 



(4.10) 
(4.11) 



where n is an integer. 

The ground and excited states for g > gc are given by 



(^-^gs,cx)'^ 

2m 



where 



-^gs.cx — ^gs,cx 



^gs , cx 



with 



^gs,ox — -J ^ 0{L^ ^), 

_ 5 , J_, .9 

n g~ gc 



0(1- 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



The case n = in Eq. ( 4.14 ) describes a delocalized 
ground state when g > gc- 

The leading term in t he en ergy eigenvalues for the ex- 
tended states (O) and ([4. 12]) is 



lim Egs cx = ^ - 
x^oo ' 2m 



(4.16) 



expl 
expj 



g/h) x] 
■g/n)\x\_ 



for X > 0, 
for X < 0, 



in the limit ij, — )■ oo. This is depicted in Fig. |[ 
Next, the excited states for g < gc are given by 



2m 



where 



(4.7) 



(4.8) 



(4.9) 



The energies of the ground state and the first excited 
state depend on g in the limit L^, — s- oo as shown in 

Fig. According to the pinning criterion given below 

Eq. ( p.20| ), the appearance of the g dependence indicates 
that these are extended states. The current carried by 
each state is thus purely imaginary. 



lim Jg, 

Lx — ^oo 
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(4.17) 



The eigenvectors corresponding to these extended states 
are 
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/ 




/ 
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/ 
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1 ^ 

, 1 , '^^c , 


^ 1 1 



Xm 

L 

(a) 







(b) 



FIG. 8. The dispersion relation Re e vs. fc((x Im e) in 
one-dimensional one-impurity case: (a) for g < gc and 
(b) for g > gc- 



+ ce 



for a; > 0, 
for a; < 



(4.18) 



in the limit ^ oo, where c is a constant. The term 
gifccxa: ^ plane wave traveling in the x direction, and 
the second term of the second is the reflection due to the 
delta potential at the origin. 

In Fig. ^, we plot the dispersion relation, that is Re e 
against fc(cx Im e). Note the level repulsion between the 
impurity level and the extended states; as we increase g, 
the localized impurity level approaches to the series of 
extended excited states, and must squeeze itself into the 
sequence. 

Finally, the solutions atthe critical point g =^ gc are 
again given by the form ([4.12 ), but with 



^gs.cx — 0(L^ InLj.), 

_ g_ 



0{L-HnL,) 



(4.19) 
(4.20) 



B. Flux-line depinning from a defect 

We now apply the above solutions to flux-line physics 
in a thin cylindrical superconducting shell with one de- 
fect in it. 

The flux line undergoes a depinning transition when 
g = gc- Its tilt slope, given as L,- — )■ cx) by the imagi- 
nary current of the ground state, vanishes for g < gc, but 
jumps to a finite value for g > gc'- 



Im J, 



gs 





g/m 



for g < 5c 
for g> gc 



(4.21) 



This directly shows that the flvix line is depinned from the 
defect at g — gc- At this depinning point, the contribu- 
tion of this flux line to the bulk transverse magnetization 
jumps by 



AM, 



Am 



(4.22) 



where A is the c ross-s ectional area, because of the rela- 
tions (|2T|) and ( |2^ ). 

Although the above jump appears to indic ate a first- 
order transition, |lj] we shall argue in Sec. VIII that 
the bulk magnetization in fact grows continuously at 
the transition when interactions are taken into account. 
There are, moreover, diverging precursors of the destruc- 
tion of the transve rse M eissner effect as 5 — > gc"- As dis- 
cussed in Subsec. II C , the normalized right-eigenvector 
gives the probability distribution of the corresp ond ing 
flux line at the top surface. The wave function (4.7) in 
the pinning regime {g < gc) defines a surface localization 
length of the flux line as £,± = (Kgs — g/h)~^, which 
diverges as g ^ gc~ , 
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9c- g 



(4.23) 



The transverse magnetization near the surface diverges 
like 

A/,, oc {x),^ (x),^/, _ ^^^^ _ ^ ^^^^ ^ 



9c -g 



(4.24) 

00), as the vortex line begins to 



(in the limit Lr,Lx 
tear free. 

We show below that the displacement of the flux line 
from the defect depends on the depth from the surface 
according to 



exp 



(4.25) 



as T — > 00, wh ere r denotes the distance below the sur- 
face. (See Sec. VIl for the many-defect case.) The pen- 
etration depth T* diverges at the depinning point g = gc 
as 



9c- g 



as g 



9c 



(4.26) 



The exponent of the penetration-depth divergence is the 
same a s the one of the surface-localization-length diver- 
gence ( 4.23 ): If we define the exponent 2; by r* ~ (C-l)^j 
we have z — 1. (Th e ex ponent is different in the many- 
defect case; see Sec. VII.) 



Th e div ergence ( 4.26 ) is derived as follows: Expand 
Eq. (2.13) with r espec t to eigenfunctions and take the 
limit — > cxD in (2.13) to obtain 



)^ ^ ^c„y" dx(V'gs I x)x{x I ^„)e-^^^"/^ (4.27) 



where 
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defect flux line 




■X 



V(x) 



bo 



-Uo 



FIG . 10. The square-well impurity potential treated in 
Sec. llVd. 



FIG. 9. The definition of the trapping length rtia 





x')dx' 


I (As 


x')dx' 



and 



(4.28) 



(4.29) 



The term ij^n — i'gs in the summation over n gives the 
position of the localization center of the ground state, 
namely the origin. Hence the leading term is due to 
the first excited state. Then the penetration depth r* 
in ( 1.25| ) is given by the inverse energy gap between the 
ground state and the first excited state, 



2mh 

9 9 

9 c ~ 9 



(4.30) 



which leads to Eq. ( [4.2(^ ). 

There are also diverging precursors on the other side 
of the transition {g ~* .gc"*") in the form of a diverging 
finite-size correction, namely the trapping length. The 
trapping length is defined as follows (see Fig. The 
average slope of the flux line for x should be given 
by the tilt slope of a flux line in the absence of defects, or 
g/m. Assume that the flux line is trapped by the defect 
at a: = over the distance rtiap- Then the total imaginary 
time necessary for a flux line to circle the periodic system 
of circumference is Lxtn/ g + Ttra,p- Thus the mean tilt 
slope (or the imaginary part of the current) should be 



given by Im Jg^ = {L^m/g + Tt 



rap J 



This, in turn, 



gives the definition of the trapping length as foUows; 



'^trap 



1 



Im Je 



m 
9 



(4.31) 



For this to be finite in the limit Lj; — > oo, the imaginary 
part of the current should have a finite-size correction of 
the form 



Im Je-. 



9_ 

m 



Ttrap.9/ m 

L'x 



(4.32) 



In fact, the correction ( 4.15| ) has precisely this form: 



Im Jtr, 



dg 

^(1 



2m 



d9 



We thereby obtain 



mh 

'^trap — ^ 



lu- 



ff 



9 \9-9c 



9- 9c 



9 -9c 



(4.33) 



(4.34) 



which has a simple pole as we approach the pinning point 
g = gc from the depinning regime. 



C. Exact solution: a square- well potential 



We can also solve the one-dimensional problem exactly 
for a single square-well impurity as shown in Fig. |lC| . 
We h ere re strict ourselves to localized states. The equa- 
tion ( A. 12 ) for the (5-potential case is now replaced by 



2kK {cosh{Lxg/'h) — cos(2fc6o) cosh[K{Lx — 26o)]) 

+ {k^ - K^) sin(2/c6o) sinh[K(L^ - 2bo)] = 0, (4.35) 



where 



^2miUa-\e\) 



and K 



V^rn\e\ 



(4.36) 



The equation for the S potential is recovered in the limit 
&o ^ and Uq ^ oo with 2boUo = Vq- 

The results are quite similar to the (5-potential case. 
The critical field of the dclocalization transition is again 
given by 



gc = lim fiK{Lx), 



(4.37) 



where k{Lx) is a solution of Eq. ( 4.35 ) with g = 0. In 
the limit ^ oo with g — Eq. ( [4.35 ) is reduced to 



2fcK cos(2fc6o) + [n^ - k^) sin(2fc6o) = 0, (4.38) 
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V. ONE-DIMENSIONAL RANDOM 
NON-HERMITIAN TIGHT-BINDING MODEL 



In this section, we present numerical results for a one- 
dimensional random system. To make the numerical cal- 
culation tractable, we use lattice non-Hermitian random 
tight-binding models (1.2), where the sites represent col- 
umn positions with random binding energies and/or hop- 
ping matrix elements, as in Fig. ^. 



A. Site-random model 



FIG. 11. The solution K„(a) of Eqs. ( [4.40D and ( [4.4lD , 
which is proportional to the dcpinning field H^c of a sin- 
gle square-well impurity. Both the ground-state function 
Ko(a) as well as the first three excited states are shown. 



or 



[Ksin(/c6o) + A:cos(fc5o)] [Kcos(fc5o) — A:sin(A:6o)] = 0. 

(4.39) 



Upon putting k = b^n and r] = bok, we are led to 



and 



where 



—2,2 -2 

K + rj — a 



K = 77 tan 77 or k = — 77 cot 77, 



boVmUo 



(4.40) 



(4.41) 



(4.42) 



is the "reduced temperature" in the sense that Ti 
corresponds to the temperature of the flux-line system. 
Ther e is always at least one solution of Eqs. ( 4.4C ) 
and (4.41), including the ground state. More localized 
excited states are available for small a, which corresponds 
to the low temperature of the flux-line system. Each state 
is localized for g < gc — fin{oo), where k is the inverse 
localization length of the state at g = 0. Figure |ll| shows 
the a-dependence of K„(a) for the ground state and the 
first few excited states. The depinning transverse field 
for the Tith state is if_Lc = (Ae / cbQ)'Kn{a) for the ground 
state as well as for the excited states. 

The divergences discussed in the previous subsection 
for the (5-function potential problem also apply to the 
present case. 



First, we discuss the site-random model. The second- 
quantized Hamiltonian in d dimensions is given by 
Eq. (|l.2| ). (Boson notation is used here, because flux lines 
behave like bosons in the delocalized regime, although 
statistics are irrelevant when the lines are strongly local- 
ized.) The hopping element is approximately pO|| 



t ~ Vbind exp(-\/2?7iVbi„da/?i), 



(5.1) 



where Vbind is a typical binding energy of the defect, and 
a is the lattice spacing. We again apply periodic bound- 
ary conditions 



for 1^ = 1,2,, 



(5.2) 



where = Li,/ a. 

The issues involved in determining the spectrum of 
Eq. are especially easy to illustrate in one dimen- 

sion. In an iV-site basis with lattice spacing a, 7i takes 
an approximately tridiagonal form, 



n 



V2 











A 



V3 



\ 



VN-1 



(5.3) 



J 



with Vn = —2Vn/t and 5 = gao/h. Delocalization of 
the eigenfunctions is traditionally associated with an ex- 
treme sensitivity to boundary conditions; the pe- 
riodic boundary conditions used here are reflected in the 
nonzero matrix elements in the upper right and lower 
left corners. For the present non-Hermitian problem, de- 
localization is reflected in an extreme sensitivity of the 
eigenvalues. Suppose the (real) eigenvalue spectrum is 
known exactly for the Hermitian site-random problem 
with g = 0. One expects all states to be localized for 
this one-dimensional problem. If the entries in the upper 
right and lower left corners are arbitrarily set to zero, it 
is easily shown that all eigenvalues for general g remain 
real and strictly equal to their values for g = 0. As will 
be shown in this section, the eigenvalue spectrum for the 
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(b) 

FIG. 12. (a) The energy spectrum of the one- 

dimensional tight-binding model with randomness A/i = 
1 and Lx = 1000a. Plots for different values of 5 = ga/h 
are offset for clarity. The same realization of the ran- 
dom potential was used for all plots here. Each eigen- 
statc is marked by a cross, (b) A blowup of a part of (a); 
Note that the real eigenvalues, corresponding to localized 
states, are independent of g. 

uniformly distributed in the range [—A, A]. For such a 
symmetric distribution, the complex spectrum is statis- 
tically symmetric with respect to the axis Re e = 0. 

Figure |l^(a) shows the d = 1 spectrum with N — 1000 
sites for various values of 5. As is discussed above, a 
complex eigenvalue indicates that a flux line occupying 



□ 



□ 

□ 
□ 



□ □ 



□ 



□ 



FIG. 14. The solid curve shows an approximate in- 
verse locahzation length determined from Eq. ( p.q ), 
plotted against Re £ for A/t = 1, = 500a and 5 = 0. 
The value was averaged over one hundred samples binned 
with energy window 0.04i. The dashed line shows an esti- 
mation of K based on the delocalization criterion = Tik. 
For this purpose, we applied increasing g to a sample of 
size Lx = 500a. The states at the mobility edges for a 
value of g {e.g. the dotted line) have the inverse localiza- 
tion length K = g/Ti. 

the delocalization criterion gc — tin, delocalized states 
appear first at the band center and the mobility edges 
move outward, in agreement with our numerical results. 

In practice, it is convenient to define the inverse local- 
ization length using the delocalization criterion gc = Tin 



|17| instead of using Eq. (5_^); if a state becomes delocal- 



ized at a certain value g^ the inverse localization length 
of the state is then k = gc/Ti, shown as a dotted line in 



Fig. |lj. The numerical estimation of k using Eq. (5.i 
difficult near the band edges, where the number of data 
points is small, and hence a large statistical error appears. 



In addition, the definition (5.8) is not tied directly to the 



asymptotic behavior exp( 

B. Probability distribution of a flux line 



Diagonalization of the lattice Hamiltonian (1.2) en- 
ables us to calculate the imaginary-time evolution of the 
wave function, and more import antly the probability dis- 



tribution of the flux line, Eq. ( 2.27 ). The expansion of 



Eq. (2.27) with respect to the energy eigenstates results 
in 

P{x\ t) ^ ^'^{tpf \ -^m) {ipni I x) X 

{x I V„) {^n I e-(^^--)-'")/''— (5.9) 



with 
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X=oo 




X=oo 



x=o 

(b) 




8>g> 



FIG. 15. The probability distribution of a flux line as 
a function of t for the one-dimensional random tight- 
binding model (1.2) with a specific random potential 
V{x) with A/t = I. The spatial size is = 100a. The 
front end of the figures corresponds to the bottom surface 
of the superconductor, while the bulk part r > 50fi/t is 
omitted, (a) Case 'g = 0.8 < with free boundary condi- 
tions, (b) Case 'g = 2.0 > with a delta-function initial 
condition. (The peak at r = in (b) has been reduced 
in size for visualization purposes.) 

2 = ^ e-^-^"/'' (V^ I Vn) (V-n I r) ■ (5.10) 
n 

In the limit L^- oo, in particular, we have 



P{x; 









^0 



n 

(5.11) 



where Ae„ = e„ - Egg. 

We used this formula to demonstrate numerically, for 
a particular realization of the random potential, that a 
flux line in the one-dimensional random system tips over 
near the surface under the influence of the transverse 
magnetic field. For g — gc, the flux line is strictly local- 
ized for large r at the strongest pinning center. As we 
increase g, kink configurations arise which allow hopping 
from one pin to the next. Figure |l^(a) shows P{x; r) for 
"free" boundary conditions at the bottom (t ~ 0) of the 
sample, corresponding to 



dx \x) 



(5.12) 



Note that the probability distribution near the surface 
has been pulled in the negative x-direction by the trans- 
verse field. Once the flux line is depinned, spiral trajecto- 
ries arise. To see this more clearly, we show in Fig. [l5|(b) 
the probability distribution for large g with the initial 
vector 



\Xo) 



(5.13) 



so that the end of the flux line is fixed to xo for t — 0. 



C. Random-hopping model 

A one-dimensional Hamiltonian with off-diagonal ran- 
domness may be more suitable for describing the exper- 
imental situation in Ref. ||l^ than the present diagonal 
randomness. Mutually parallel twin boundaries might 
represent the pinning potential in this situation. As dis- 
cussed above, we project out the coordinate which is par- 
allel to the twin boundaries and perpendicular to the r 
axis. Randomness arises from the separation of the bind- 
ing twin boundaries rather than from the strength of the 
pinning. Thus we have 



(5.14) 



where tj = Vbind exp[(— A ± g)aj/?i], A = V2wT4ind, and 
ttj is the separation between the jth and {j + l)th binding 
impurities. If the twin boundaries are located randomly, 
the random separation follows the Poisson distribution 
P{aj) = a~^e~°j/'', where a is the average separation. 

The delocalization phenomenon discussed above arises 
in this off-diagonal case as well; see Fig. |l^. In this ex- 
ample, we neglect the randomness embodied in the factor 
exp[±gaj/?i] and use a square distribution of the hopping 
elements: 

where i is a random variable with probability 
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FIG. 17. The energy spectrum of the two-dimensional 
non-Hermitian tight-binding model without impurities, 
with gx — 9y ~ 1-0 X h/a and = Ly ~ 20a. 



27r 47r 
fc, = 0,± — ,± — , 



(6.3) 



or 



Re e 



cosh(ga/7i) 



Im e 



n 2 



sinh(ga/?i) 



cos^ 



(6.4) 



Thus the spectrum consists of ellipses with various radii 
as is shown for = Ny = N = 20 in Fig. |l^. Two levels 
with and ky interchanged are degenerate at each cross 
of the figure except for the point e = 0, where N levels 



with kx 



7r/a (mod 27r/a) are degenerate. The 



eigenf unctions have the usual Bloch form: 



Ifc) =6tlO) = 



1 



^e*-6t|o). 



(6.5) 



The important phenomenon of level repulsion in the 
complex plane can be illustrated with one attractive 
point impurity: 



Vx = -Uo5xxo, (6.6) 

where Xq denotes the position of the impurity. Figure [l8| 
shows the results for the 20 x 20 lattice with periodic 
boundary conditions. Without the transverse field, we 
have one impurity level separated from a cluster of (de- 
localized) excited states, as is seen in Fig. |l^(a). As we 
introduce g diagonally {gx = gy), most excited states 
acquire imaginary eigenvalues (Fig. |l8|(b)), showing a 
spectrum similar to the impurity-free case. Fig. QTl The 



FIG. 18. The energy spectra of the two-dimensional 
non-Hermitian tight-binding model with one attractive 
impurity of well depth Uq. We show here the case 
Lx — Ly = 20a and C/q = lOi with (a) Qx — gy ~ 0, 
(b) gx ~ 9y = 1-0 X h/a, and (c) gx — gy — 2 x fi/a. 



localized impurity level does not change its energy as long 
Bs\g\ < gc, while the extended excited levels expand, fol- 
lowing Eg . (3.4) of the impurity-free case approximately. 
In Fig. [l8|(c), we observe level repulsion between the im- 
purity level and the excited levels. Note that the two- fold 
degeneracy of the pure system is split on the side of the 
spectrum, where the impurity state is about to enter. For 
even larger the bound state enters the region of de- 
localized levels and the spectrum is again elliptical with 
an extra state near the origin. 



B. Random case 



In the two-dimensional random tight-binding model, 
we again find energy bands of localized levels bounded 
by a mobility edge. For certain values of g, however, ex- 
tended and localized states are mixed in a complicated 
way near the band center. It appears that we have 
three regimes with respect to the field g. One is the 
pinning phase for small g. Another is the Bloch-wave 
regime for large g. Finally, there is an exotic intermedi- 
ate regime with chaotic eigenvalue spectra. We describe 
these regimes on the basis of our numerical data. 

First, it is widely accepted for g = that all eigen- 
statcs of two-dimensional random systems are localized 
with finite localization lengths (except for a possible ex- 
ception at the band center |^). Hence there is again a 
region of small g where all states remain localized. As 
discussed earlier, these localized states retain the real 
eigenvalues which they had for g = 0; see the energy 
spectrum in Fig. |l^(a) for an example with Nx = Ny — 
40. 

As we increase g, delocalized levels with complex eigen- 
values appear as is shown in the energy spectrum in 
Fig. |l^(b). A remarkable difference of the energy spec- 
trum from the one-dimensional case is that localized 
levels and delocalized levels coexist in the same range 
of Re £ in the two-dimensional case. This phenomenon 
can be understood in terms of anisotropic localization of 
the Hermitian system (in the case g ^ Q). To take ad- 
vantage of anisotropic fluctuations in the site potential, 
there can be states with different localization lengths in 
specific directions. As shown schematically in Fig. in 
a large system there can be another state with nearly the 
same energy, but whose contours are rotated by 90° from 
the first one. According to the argument developed in 



Sec. II], the delocalization point of each state is gc = 

(c) 
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(d) 

FIG. 19. The energy spectra of the two-dimensional 
non-Hermitian tight-binding model with site random- 
ness. We show here the case ~ Ly = 40a and A = lOt 
with (a) g-x ~ Qy — 0.7 x fi/a, (b) gx = gy = 1-2 x h/a, 
(c) gx = gy ^ 1-5 X h/a, and (d) = gy ^ 4.0 x h/a. 



where k is now the inverse localization length of the state 
in the direction of g. The state with the largest local- 
ization length in the direction of g gets delocalized first. 
Up to the delocalization point, however, both states have 
nearly identical real energies. 

After passing through an intermediate region whose 
energy spectrum is exemplified in Fig. |l^(c), we move 
onto the region of larg e g. The energy-spectrum struc- 
ture shown in Fig. |l9| (d) is similar to the impurity- free 
case shown in Fig. |17| . 

It is tempting to conclude from Fig. |l^(d) that the 
large-g limit is well described by the extended Bloch 
wave functions, just as in one dimension. However, unlike 
d = 1, the Bloch approximation breaks down in d = 2 
for sufficiently large systems for any value of g. To see 
this, recall that all Bloch states \kx, ky) with kx ^ ky are 
two-fold degenerate in the impurity-free case: If g points 



along the lattice diagonal as in Eq. (6.2), then jfe^;, fcj,) and 
\ky,kx) have the same complex energy. More generally, 
states related by a reflection across the g axis are degen- 
erate. (In d = 3, the degenerate states lie on a circle in 
{kxTky,kz)-^]i&ce centered on a line passing through g.) 
This degeneracy is split by the randomness. Degenerate 
perturbation theory requires diagonalization of the set of 
2x2 matrices, 



£o{kxj ky) V{kx ky, ky kx) 

V{ky kx, kx ky) EQ^kx, ky) 



(6.7) 



where 
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FIG. 20. A schematic view of two localized wave func- 
tions of the Hermitian system in two dimensions. The 
sohd curves indicate contours on which the locahzed wave 
function takes the same value. The anisotropy arises be- 
cause of anisotropic fluctuations in the impurity distri- 
bution. We show an impurity fluctuation which leads to 
two nearly degenerate wave functions related by a 90° 
rotation. 



Ijrp J_J, 



(6.8) 



x,y 



is the Fourier-transformed disorder potential and 
£o{kx, ky) is given by Eq. ( |6.2[ ). The eigenvalues of 
Eq. (@ are 



where 



£±{kx,ky) £o{ka;,ky) ± A, 



(6.9) 



The real splitting of the degenerate doublets predicted 
by Eq. (3^) is clearly visible in Fig. [ig|(d). 

The Bloch spectrum is only a good approximation pro- 
vided these corrections are small compared to the spacing 
between the doublets. However, a typical spacing be- 
tween the real pa rts of the Bloch levels for ~ Ly ^ L 
given by Eq. (6.2) is 



Ae 



dRc So 



dk 



Ak : 



2tacosh[^j—. (6.11) 



The splitting A due to the randomness is thus only a 
small correction to the Bloch levels whenever 



47r 



ta cosh 



(f) 



> A, 



(6.12) 



Fig. |l^(c). The Bloch states are nondegenerate in d = 1, 
so this problem does not arise. 

Some insight into the meaning of the chaotic eigenvalue 
spectra follows from tracking the ground-state wave func- 
tion as a function of g. As shown in Fig. the ground 
state first streaks out in the direction of g. (For very 
large g, the wave function eventually broadens to cover 
the entire lattice. However, as discussed above, we be- 
lieve that this Bloch-wave behavior is an artifact of the 
small system size.) Arguments given in Ref. ^ show 
that the ground state can never simply delocalizc in the 
g direction while remaining localized in the perpendic- 
ular direction when L oo. In Ref. 0, it is argued 
that the long-wavelength, low-frequency behavior of the 
non-Hermitian Schrodinger equation in 2 + 1 dimensions 
for large g is described by a 1 -I- 1 dimensional Burgers' 
equation. The tilted flux lines described in this way wan- 
der away from the g direction to take advantage of ex- 
ceptionally deep minima in the disorder potential. This 
mapping predicts a streaked-out but roughened ground- 
state wave functions with roughness exponent 2/3. Q At 
present, our system sizes are too small to provide a good 
quantitative check of this hypothesis. 

For system small enough so that Bloch states are a 
good first approximation for large g, most states are 
approximately two-fold degenerate as discussed above. 
However, for an L x L square lattice, there is a very 
large L-fold degeneracy for the state at the origin of the 
complex energy plane. Perturbation theory explains the 
removal of the degeneracy in the following way. We show 
that the first-order perturbation theory reduces to a di- 
agonal one-dimensional random model, whic h is readily 
solved. We take the hopping term of Eq. (1.2) as the 
non-perturbative part and the random-potential term as 
the perturbation so that the zeroth-order spectrum may 
be given by Eq. (6.2), or Fig. |l^. The zeroth-order wave 



(6.10) 

functions \kx, ky) of the degenerate levels in question sat- 



isfy the relation kx — ky — tt / a (mod 2n/a). The secular 
equation for the first-order perturbation of the degener- 
ate levels consists of the matrix elements 



kx J ky 



Y.^Abx 



^XJ 



(6.13) 



x,y 



where we used ky — ky = kx — k'^ (mod 27r/a). The above 
matrix elements are, in fact, equivalent to the momentum 
representation of the one-dimensional Hamiltonian 



(6.14) 



a condition which is always violated in suffi cient ly large 
systems as L — > oo. When the inequality ( 6.12| ) is not 
satisfied, level repulsion interacts with the randomness 
in a complex way to produce chaotic spectra like that in 



where 



(6.15) 
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(d) 

J _ FIG. 21. The ground-state right-eigenfunctions in two 
I dimensions. We show the case = Ly = 40a and 
A — lot with (a) gx — gy — 0, (b) gx = gy = 1-7 x h/a, 

(c) gx = 9y = 3.0 X h/a, and (d) gx — gy — 5.0 x h/a. 
The vertical scale was rescaled by 1/10 and 10 in (b) and 

(d) , respectively. The highest peak in (b) was reduced in 
size for the visuahzation purposes.) 



X 



3.0 



with = X + y (m od L ) and rj = x — y (mod L). The 
effective potential ( |6.15 ) is thus obtained by integrating 
the original potential Vx^y along a line parallel to g. The 
first-order perturbation spectrum is readily solved, since 
the effective Hamiltonian ( |6.14|) is already diagonalized. 
Thus we arrive at the first-perturbation energy 



-(1) 



(6.16) 



for ^ = 0, 1, . . . , L — 1, which accounts for the horizontal 
spread of eigenvalues in Fig. f9|( 



VII. TRANSVERSE MEISSNER EFFECT AND 
ITS PENETRATION DEPTH 



X 



In this section, we discuss the penetration depth asso- 
ciated with the transverse Meissner effect for supercon- 
ductors with columnar defects in the localized regime. 
The vortex density is assumed to be finite. The surface 
deflection which defines the penetration depth for a single 
vortex line is shown in Fig. |l^(a). We discuss only field 
penetration in directions parallel to the columns. For 
a more general discussion of various penetration depths, 
frequency dependence s, etc. , see Ref. |23] . 



As shown in Subsec. IV B for one-dimensional case, the 



displacement of the pinned flux line near the surface of 
the superconductor with one defect has, in general, an ex- 
ponential r dependence. The penetration depth r* has 
the singularity T* ~ ^ (5^5c)^^- We show below that 
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the mean displacement at depth r in the many-defect 
case in d dimensions is given by a stretched exponential 
form 



(((^ 



exp[ 



^a{r/T*)'/'-''+'^], (7.1) 



where we now assume for simplicity that the field g is 
parallel to the x axis. The average displacement in direc- 
tions perpendicular to g vanishes because of statistical 
symmetry. The quantity (x) oo is the center of the local- 
ized state in the bulk, which is identical to the average 
position of the corresponding state for g — 0. Moreover, 
we show that the penetration depth has the singularity 
T* ^ ^ (gc—g)^^ with the d yna mical exponent z = d, 
thus justifying Eqs. (1^) and (|l.9|). Note that d is the 
dimensionality of the quantum system corresponding to 
flexible lines in d -I- 1 dimensions. 

Consider a low but finite density of interacting flux 
lines. As discussed in Sec. ||, we fill up the localized 
states in order of increasing energy up to the average 
chemical potential e = ji. Consider the deflection of the 
most unstable (pinned) flux line near the surface. Since 
we forbid double occupancy of localized states, we can 
approximate 



c„ / d''x{i^^ I x)x{x I V„)e-"^^"/^ (7.2) 



where 



Ji^f, I x')d'ix'' 



and 



(7.3) 



(7.4) 



in contrast to Eqs. ( ^.27D , ( [4.28| ), and ( [4.29[ ) appropriate 
to a single line. Here ipf^ denotes an eigenstate at the 
chemical potential, and the summation is restricted 
to states with the energies e > /x. Note that the coeffi- 
cient Cn for an extended state will be quite small owing 
to oscillatory factors in the integration over the space. 
Because of the exponential factor in (7.2), the main con- 
tributions come from localized states near and above the 
chemical potential. We use the as ymp totic form ( |3.5| ) 
for localized states to estimate Eq. (7.2). The term with 
"0" = V'm gives {x)^. The other leading terms after the 
integration are approximately 



Mr) 



(7.5) 



where 




FIG. 22. A kink in a vortex configuration near the top 
surface of a superconducting sample. 

r„ = \xn - Xf,\, cosOn = 9 • {Xf, - Xn)/{\g\rn), and 
is the summation over localized excited states with 
the energies e > fj,. The inverse localization length 
was approximated by k(/j.). The new coefficient c„ is 
Cn ~ c„r„ cos^„. As in conventional Mott variable-range 
hopping in semiconductors, we estimate the energy 
difference Ae„ from [[lO[ 



5(M)r^Ae, 



1, 



(7.7) 



where g{fJ-) is the density of states at the chemical poten- 
tial. 

The quantity /n(T) gives the energy of a tilted flux- line 
configuration; see Fig. The first term is the energy 
due to a kink joining columns at x^ and a;„. Hence 
the energy is proportional to the width r„ of the kink. 
The second term is the energy loss arising because the 
flux line stays at Xn over the distance r rather than at 
the most stable position a;^ . T his energy loss is propor- 



tional to T. Because of Eq. (7.7), the first and the second 
terms compete: The further the flux line hops from x^, 
the more the kink energy costs, but the lower the bind- 
ing energy at Xn- In variable-range hopping of electrons 
in semiconductors, Ae„ is the energy difference be- 
tween localized electronic states, and the kink energy is 
replaced by a WKB tunnelling matrix element. 

Figure ^ shows the relaxation function (x)^— (x)^ for 
the one-dimensional tight-binding model with a particu- 
lar realization of randomness. In this example, only two 
state s co ntribute significantly to the summation in 
Eq. ( [7.5| ) . Upon averaging over many realizations of ran- 
domness (or close to the depinning transition), we expect 
contributions from many states as depicted schematically 
in Fig. [23. For a fixed r, the largest contribution to the 



summation in (7.5) comes from the state at 



rd 



hg{fi)[nK{fj.) - g] 



i/(d+i) 



(7.8) 



/„(r) = {hK{n) - gcos9n)rn + rAen/Ti, 



(7.6) which yields 
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FIG. 23. The r dependence of the displacement of 
a flux line for a realization of randomness in the one- 
dimensional tight-binding model. Since only two terms of 
/„ contribute to the summation in ( [7.5| ) in this exam- 
ple, the whole curve approximately consists of two lines. 
The slope of each line is — Ae„/?i, while its y-intercept is 
— {hKn — g)rn- The parameters are Lx — 500a, A/t = 1, 
and g/t — 0.6. 



mm Jn = a\ — , 



with 



{h>,{fi)-g)d- 



(7.9) 



(7.10) 



and a — 4- 1/d). This gives the stretched ex- 



ponential form (7.1). Since the displacement at the sur- 
face, or the surface transverse magnetization, is given by 
^± ^ M±s ^ i^K — 5)""^, we have r* ~ 

We obtain the stretched exponential form only in the 
limit T 00. Deep in the Bose glass phase, it may 
be difficult to observe this form experimentally, because 
the amplitude of the surface displacement is small in this 
limit. However, the characteristic scale r* in the r di- 
rection diverges, r* ^ {'hn{^) — g)^"^, where stretched- 
exponential relaxation should be relatively easy to ob- 
serve. 



VIII. EFFECT OF INTERACTIONS IN THE 
DELOCALIZED REGIME 

A. 1 + 1 dimensions: Continuous phase transition 

We conclude with a discussion of interaction effects. 
The delocalization transitions discussed above were 
treated using an independent particle (or "independent 
vortex line") picture. As outlined above, interactions 



FIG. 24. Schematic of the situation where many dif- 
ferent exponential terms contribute to the summation 
in (IJ). 

can be taken into account in the Bose-glass phase, pro- 
vided that we forbid multiple occupancy of the localized 
states. The physics here is similar to the bands of lo- 
calized impurity states describing electrons in disordered 
semiconductors. |2^] In both Bose and Fermi glasses, lo- 
calized states are filled in order of increasing energies 
up to a chemical potential. In the localized regime, the 
differences between Bose statistics of repelling flux lines 
interacting with columnar defects and Fermi statistics of 
electrons in disordered semiconductors are not expected 
to be important, jl^ 

Interactions must be handled differently in the delo- 
calized phase. Consider what happens to the states de- 
scribed by the (1 -I- l)-dimensional non-Hermitian spec- 
tra of Fig. [ij (a) with increasing field iJy along the r 
axis. We assume that the tilt field lies in a range such 
that there is a mobility edge separating low-energy local- 
ized states from high-energy delocalized ones. The field 
-ff II controls the chemical potential of the equivalent dis- 
ordered boson system. [gO| As i/|| increases, we fill the 
unoccupied levels in order of increasing energy to obtain 
the ground state. Eventually all states below the mobility 
edge are filled, and additional vortices must then go into 
extended states above this boundary. These extended 
states describe macroscopically tilted vortex lines, and 
as we can see from Fig. ^ the corresponding tilt slope 
is finite at the mobility edge. Interactions have a weaker 
effect on delocalized, tilted lines, so we expect that many 
lines can be accommodated by only a few delocalized 
states with energies just above the mobility edge. In the 
presence of thermal fluctuations, there is not a sharp dis- 
tinction between tilted and untilted lines; any individual 
line will be both tilted or localized in different regions 
along the r axis. Alternatively, we can imagine that the 
delocalization transition happens at fixed TJy , with in- 
creasing external tilt field H± , as in Fig. ||. 

A physical picture of the delocalized phase in 1 4- 1 
dimensions in the presence of interactions has been de- 
veloped by Hwa et al. pS]. Tilted lines are represented 
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by chains of kinks, and the density of these chains goes 
continuously to zero with decreasing tilt field. Thus we 
predict that a continuous phase transition is possible, 
in contrast to arguments based on depinning of a single 
line, which suggest a first-order transition, jl^ A phe- 
nomenon like Bose condensation describes the physics of 
the tilted fraction of the lines, consistent with many lines 
entering just a few extended states above the mobility 
edge. However, because phase fluctuations are so strong 
in 1 -|- 1 dimensions, correlations of the boson parameter 
decay algebraically to zero at large distances instead of 
approaching a nonzero constant. There is also algebraic 
order in a translational order parameter associated with 
the density of tilted lines; the tilted phase is in fact a 
"supersolid." E^] 



B. 2 + 1 dimensions: Existence of a phase transition 

The physical picture of lines in the delocalized regime 
is more complicated in 2 + 1 dimensions, as suggested 
by the intricate spectra for finite g displayed in Fig. pj| . 
However, here again there is a range of intermediate tilt 
fields where only localized states exist at low energies. 
A "mobility edge" now separates real eigenvalues in the 
localized regime from a region along the real energy axis 
where extended and localized states coexist. Building 
up a ground state as in 1 -I- 1 dimensions, we would ex- 
pect some vortices to enter extended states describing 
tilted lines once the chemical potential (controlled by 
i?l|) crosses this mobility edge. A bosonic "superfluid" 
fraction of tilted lines will coexist with a "normal fluid" 
fraction of lines localized on columns. As in 1 1 dimen- 
sions, however, any single line will participate in both 
"fractions" as it crosses a sufficiently thick sample. Al- 
though a single tilted line probably goes into a glassy 
ground state described by Burgers' equation, inter- 
actions in a dense liquid of tilted lines can screen out 
random pinning potentials in a weakly perturbed "su- 
perfluid" phase of entangled lines. [^. The tilted liquid 
in this dense regime should be pinned only weakly, and 
exhibit a linear resistivity. 

If the tilted regime is entered by increasing H± at fixed 
one is shifting the "mobility edges" in the (2 -I- 1)- 
dimensional spectra at fixed chemical potential. The se- 
quence of possible phases probably resembles vortex mat- 
ter subjected to point disorder near Hd. |26| For small 
H± , the Bose-glass with its transverse Meissner effect re- 
sembles the usual Meissner phase ioi H < Hd- Suppose 



now that H± is increased until H± > H±c, such that 
the lower mobility edge drops below the chemical poten- 
tial. As illustrated in Fig. ^, one probably encounters a 
flux-liquid phase (with a finite fraction of tilted, entan- 
gled lines), if H± is increased at high temperatures. At 
low temperatures and large tilt fields, however, the sta- 
ble phase is probably a crystalline vortex-smectic phase, 
where the tilted lines are arranged in sheets periodically 
spaced along the column direction to minimize the in- 
teraction energies. Hwa et al. |^ have discussed a 
theory of how a continuous transition directly from this 
vortex smectic into the Bose glass might proceed. How- 
ever, it is hard to rule out the possibility of a melted 
silver of tilted flux liquid at all nonzero temperatures 
when the density of tilted lines becomes very small near 
the transition curve H±ciT), similar to what happens in 
vortices subjected to point disorder near Hci{T). It is 
not yet known if a point disorder is sufficient to produce 
a distinct "vortex glass" phase with a finite vortex den- 
sity interposed between the flux liquid and the Meissner 
phase. |3| There are similar uncertainties about a dis- 
tinct glassy phase near H±ciT) in the present problem. 

Although we will not resolve these uncertainties here, 
we can show that there must be at least one phase tran- 
sition separating a Bose glass at small H± from a tilted 
"superfluid" of entangled lines when H± is large. In 
other words, the transition suggested by the existence of 
a sharp "mobility edge" in the (2-1- l)-dimensional single- 
line spectra survives the imposition of repulsive inter-line 
interactions. We proceed by first noting that the tilt 
modulus, 



C44 



47r \dH\ 



(8.1) 



is divergent in the Bose-glass phase which results for 
small H± . [0 This infinity is a consequence of the trans- 
verse Meissner effect: Although the lines tilt near the 
sample surfaces in response to a perpendicular external 
field (as in Fig. |2^ ), there is no macroscopic response in 
an infinite sample, so dB±/dH± = 0. We will show that 
the tilt modulus C44 remains finite in a "superfluid" liq- 
uid of entangled, tilted lines which results for H± > H±c, 
so that there must be at least one sharp phase transition 
as one increases H±. The tilt modulus is proportional 
to the inverse superfluid density of the equivalent boson 
system, and the corresponding results for the superfluid 
density ps are indicated in Fig. ||. 



We start with the generalization of Eq. (2.1) for N 
lines {xj{T)}, namely, ||l^ 



Efiux [{xj{t)}] = yJ2 



dT 



dx., 



dT 



dxi 



dT 



/ ^ dTV^,{\x,{T)-X,{T)\), (8.2) 
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where g = (j>QH j_/ {Air) and we have added a repulsive 
potential Vint between the lines. We neglect point disor- 
der and assume for simplicity that the columnar pinning 
potential is weak or the temperature is fairly high, so 
that the transition out of the Bose glass occurs for rela- 
tively small values of g cx H± . The assumptions of small 
line tilts and local (in r) interactions between flux lines 
which justify Eq. (^^) [|lO| are then satisfied even for 
H±^ > Hj^ c- T he many-line partition function associated 
with Eq. (O) is 



N 

2(9)^11 J -Dx.ir) 



--Etlux[{a!j(T)}]/(fcBT) 



(8.3) 



Upon defining an operator Xop = Y^^=i ^'^^i'^S 
position eigenstates |a;i, . . . , a;jv) such that 

jN \ 

Xop \xi,..., xn) = \^Xj \ \xi, . . . , xn) , (8.4) 

standard manipulations allow us to write this multidi- 
mensional path integral in terms of a quantum mechani- 
cal matrix element, 



= I cPr\ ■ ■ ■ d?r'j^ / (fn ■ ■ ■ (frN (r[--- dr'j 



N 



gg-x,^/{kBT)^- 



H{0)LJ{kBT) -g-Xa^/{kBT) 



ri ■ ■ ■ dr 



N 



(8.5) 



where the "Hamiltonian" 7i(0) is 
{ksT) 



2ei 



int ( l-^i ^ ■^j I ) ■ 



(8.6) 



Upon noting that 



gg-Xop/(fcBT)g--H(0)L,/(feBT)g-g-Xop/(fcBT) 



- V —p9-^°p/i''BT) . ____ . 

71 — ^ 



= ^-H(g)L^/{kBT) 



(8.7) 



where 



Hig) = e9-^-/('=«^)H(0)e-s-^op/(feBT) 

-I / . rp \ 2 N 



int \\Xi Xj 



Upon referring to Table |, we see that this Hamiltonian 
is the generalization of Eq. (^]^) for many interacting 
vortex lines. 

Note that if a many-body eigenfunction 
^'„(a;i, . . . , a;jv; 0) of 7i(gr = 0) for localized lines in 
the Bose glass phase is known, a potentially exact right- 
eigenfunction with the same energy for g ^ is then 

{xi, ...,XN;g)^ ^g-x.^/ikBT)^^ (a^i, . . . , xn; 0) 

(8.9) 

as follows immediately from the first line of Eq. (^.^). 
As in the single-particle case, this gauge transformation 
connecting the g = and g problems only works 
provided that the new eigenfunction is normalizable. 

To treat the tilted phase of entangled lines, it is conve- 
nient to modify the original path-integral partition func- 
tion via the change of variables. 



Xjir) = x'Jt) + ^T. 



£i 



(8.10) 



Upon transforming to a quantum-mechanical matrix el- 
ement as before, we arrive at 



z{g) 



■d^r' 



N 



d^n 



■d^rN 



N 



T exp - 



n'{g;T)dT 



ri---rN 



(8.11) 
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with 



H'(g;r) 



2£i 



N 



N 



.12) 



where the constant imaginary vector potential is missing, 
but V[xj] has been replaced by a r-dependent disorder 
potential 



V' [xj,t] = V 



£l 



(8.13) 



The symbol T in front of the exponential in Eq. ( 8.11 ) 
stands for time ordering. The above new Hamiltonian 
(related to Eq. ( |8.8| ) via an imaginary Galilean transfor- 
mation 1 13 1) describes a set of vortex lines moving along 
the T axis in the presence of a set of parallel, tilted colum- 
nar defects. The response functions for a liquid of inter- 
acting entangled lines in the presence of tilted disorder 
have been discussed and reviewed by Tauber and Nelson. 
pof The vortex tilt moduli both parallel and perpendicu- 
lar to the plane of tilt in this case are different. However, 
both are explicitly found to be finite, in contrast to the 
infinite tilt moduli in the Bose glass. Thus there must 
indeed be at least one genuine phase transition for the 
(2 -|- l)-dimensional system with increasing H ± cx g in 
the presence of interactions and disorder. 
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APPENDIX A: EXACT SOLUTION OF THE 
ONE-DIMENSIONAL SCHRODINGER 
EQUATION WITH A POINT IMPURITY 

In this Appendix, we describe the derivation of the 



single-impurity results given in Subsec. W A . The deriva- 
tion is subtle even in the Hermitian case, since we are in- 
terested in finite-size effects as well as the thermodynamic 
limit. Throughout this Appendix, we assume g > with- 
out loss of generality. 

We have found it useful to solve the Schrodinger equa- 
tion Hip^{x) — eip-^'{x) in a way different from that 
sketched in Ref. H. Because of the periodic boundary 



conditions ( |4.2| ), we define the Fourier transformation in 
the form 



(A.l) 



where the summation runs over k = 2Tm/Lx for integer 
n. The inverse transformation is given by 



-'^""dx. (A.2) 
The Schrodinger equation becomes 

-^^■(0) = 2meij'\k). (A.3) 



{hk + i5)2^^(fc) - _ o^.J.flv 



Upon assuming that C = "(/"^(O) is finite, we have 

-1 



v;^(fc) 



mVbC 



(A.4) 



where K{e) = \/2m£/fL. The real-space wave function is 
given by substituting 'ij)^{k) in Eq. ( [A.l| ) with Eq. (A.4). 

The condition ■i/;^(0) = C results in the equation for 
the energy spectrum e; namely 



2tt 



k + z^ 
n 



K\s) 



(A.5) 



Let us first analyze this equation in the limit Lx ^ oo: 



dk 



k + i- 



mVo 



(A.6) 



The integrand has two poles at fc = ±K — ig/h in the 
complex plane of k. For |Im K\ < g/h, the integral ( |A.6| ) 
vanishes when we close the integration contour in the 
upper half plane. Thus, nonzero solutions arise only for 
|Im K\ > g/h. In the following, we derive solutions for 
|Im K\ > g/h and for |Im K\ = g/h separately. 

For |Im K\ > g/h, the poles straddle the real axis, and 
we have 



TTl 



TTh^ 

mVo' 



or 



K 



.mVo 



(A.7) 



(A., 



Thus we arrive at the unique localized ground-state solu- 
tion for |Im K\ > g/h given as Eq. (4.4). The boundary 
for the existence of this localized state is the critical field, 
which leads to Eq. (^). 

For Im K — g/h, one of the poles of the integrand of 
Eq. (A.6) approaches the real k axis, which makes the 
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FIG. 2 5. Th e integration contour used in the evaluation 
of Eq. (A.IC). The crosses on the horizontal axis indicate 
the poles of the factor (e*^^ — while the other two 

crosses indicate the poles of the function R{z). 

evaluation of the integral in the limit ^ oo difficult. 
Hence we return to Eq. (A. 5), and explore the solutions 
including finite-size corrections, using the formula 



2lT 

T 



E„ / 27r 
R f —n 



n— — oo 



L 



^Res 

C 



R{z) 



2m 



-;C 



(A.9) 



where R{z) is a rational function with the conditions that 
(i) the order in z of the denominator is greater than the 
order of the numerator at least by two, and (ii) the func- 
tion does not have any poles at integral points. The 
summation with respect to C, runs over the residues aris- 
ing from all poles of R{z). The formula follows from 
evaluating 



dzR{z) 



27ri 



z,iLz 



- 1' 



(A.IO) 



with the integration contour as in Fig. Thus Eq. ( |A.5| ) 
becomes 



iKL^+L^g/h _ 2^ 



mVo 
(A.ll) 



K 



cosh (^La-^ — COS (LxK) 



sin {L.,K) - 0. 



(A.12) 



The wave function corresponding to each solution of 
Eq. (|A.11|) is 



^iKx-\-xg jh 



^ — iKx-\-xg /h 



iKLj,+Lj:g/h _ I ^-iKL^-\-Lj,g/h _ ^ 



(A.13) 



for < X < Lx- In the case of limi^_»oo Im K{Lx) > f , 
putting K to be pure imaginary results in the ground- 
state solution (A.8) and ( fl.Tj ) for g < gc in the limit 
Lx — > oo. 

We are now in a position to discuss the case 



lim ImKiLx) = |, 



(A.14) 



where K{Lx) denotes a solution of Eq. ( A.12 ). This case 
includes all delocalized wave functions. By setting 



with 



lim 5k,{Lx) = 0, 



Eq. (A.ll) becomes 



T 9 - 
n 



- (f - f ) 



in the limit ^ oo, where 

w = Lx6k{Lx) — iLxk[Lx). 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



We notice that n{Lx) and k{Lx) have the same order of 
magnitude with respect to L^. 

Fo r g ^ gc, the second term of the right-hand side of 
Eq. (A.17) is negligible compared to its first term, and 
hence 



lu- 



ff 



9 -9c 



(A.19) 



This g ives a series of solutions in addition to the ground 
state (A^)- For g < gc, we have 



Lx6k = In , 

gc - 9 

Lxk = In(-l) = i7r(2n + 1) 



(A.20) 
(A.21) 



with n a n inte ger, i.e. the excited states given in 
Eqs. ( 4.10 ) and ( 4.11 ). For g > gc, we have 



9 



LxSk = In ■ 

9 - 9c 

Lxk = In 1 = 2z7rn, 



(A.22) 
(A.23) 



or the excited states given in Eqs. ( |4.14| ) and ( |4.15| ). 

Finally, for g = gc, E g. ( lA.lTf) yields a solution w = 
0{\nLx), or Eqs. (Ol and ([4.20D. 
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